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Due to the time reversal invariance of the angular momentum operator J^, the average energies 
and variances at fixed J for random two-body Hamiltonians exhibit odd-even-J staggering, that 
may be especially strong for J — 0. It is shown that upon ensemble averaging over random runs, 
this behaviour is reflected in the yrast states. Displaced (attractive) random ensembles lead to 
rotational spectra with strongly enhanced BE2 transitions for a certain class of model spaces. It is 
explained how to generalize these results to other forms of collectivity. 
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Even-even nuclei are more bound than their odd neigh- 
bours (OES, odd-even staggering). Their ground states 
have always spin J = (JOD, J = dominance), and 
with excedingly few exceptions their first excited state 
has J = 2 and decays through an enhanced BE2 transi- 
tion. These systematic features have been traditionally 
interpreted in terms of specific components of the nu- 
clear interaction. OES, for instance was attributed to 
the pairing force but it has been argued that, in small 
Fermi systems, it is a universal phenomenon associated 
with deformation Though the claim is too strong, 
there is some truth in it as it can be viewed as a con- 
sequence of the Jahn- Teller effect which produces OES 
through the filling of time-reversed pairs in doubly de- 
generate single particle states, e. g., Nilsson orbits 

The discovery that random interactions lead to J = 
ground and J — 2 first excited states far more often than 
statistically expected ^ has brought the subject of sys- 
tematic features into sharper focus. In particular, OES is 
often present, even in the absence of either pairing forces 
or deformed fields , which means that it should not be 
associated to a specific form of coherence. We are left to 
find general causes. 

In what follows it will be argued that time reversal (T) 
invariance is at the origin of OES and JOD. Later, it will 
be shown that the attractive nature of the interaction 
leads to BE2 coherence so far detected in a randomized 
IBM context ||^ but not in shell model (SM) simulations. 

Before going into the simulations themselves, we define 
notations and collect some basic results and predictions 
based on T invariance. 

We shall rely for insight on the low moments of the 
Hamiltonian Hat fixed J, which determine the level den- 
sities Calling dj the dimensionality of the vector 
space, the centroid is Ecj = dj^ti{H)j = dj^{H)j, and 
the variance a"^ = dj^ {{H — EcjY)j- 

Using r, s, . . . for subshells and i = rr^ ■ ■ ■ for individual 
orbits, H can be written in uncoupled or coupled form, 
in terms of pair creation and annihilation operators 

i<j,k<l r<s,t<u,J 



The W;4tti matrix elements will be taken to belong to 
the Gaussian orthogonal ensemble (GOE), i.e., to be real 
and normally distributed with mean zero and variance v'^ 
for the off-diagonals and 2v^ for the diagonals. For n > 2 
the Hamiltonian matrices are said to belong to the two 
body random ensemble (TBRE). An alternative should 
be considered: GUE, where U stands for unitary, in which 
case the matrix elements are complex (general Hermitian, 
rather than symmetric, matrices). TBRE derived from 
GOE Hamiltonians are T-invariant, while those derived 
from GUE ones are not ||^. 

In Ref. iQ it was found that both ensembles lead to the 
same behaviour for OES and JOD. At first sight the result 
is surprising since T invariance is at the origin of such 
behavior, as we shall see soon. The paradox is resolved 
by noting that the operator is T invariant, and states 
of good J, transform as T\JM{^ (-)'^+^|J - M) [|. 
At M — 0, which contains the information on all states, 
eigenstates of are also eigenstates of T, with eigen- 
value {—)"''■ true for TBRE derived from either GOE or 
GUE. Now examine the consequences. 

In a determinantal basis, at M = 0, if a state and its 
time reversed are different, the resulting even and odd 
combinations (e.g., \jm j' — m) dz^—y'^^ \j'm j — m)) con- 
tribute equally to the traces of H''. The self conjugate 
states are always even (e.g., \jmj — m)). Therefore, at 
fixed number of particles n there are always more even J 
states than odd ones. 

Obviously, the difference in average energies between 
even and odd states, Cj will depend on the sign for the 
average over selconjugate states. Upon ensemble aver- 
aging ^ will vanish, and we expect no staggering between 
Ecj even and odd states. It is equally obvious that for 
the variances all contributions are positive. Hence the 
average aj will always be larger for even states than for 
odd ones, and staggering will persist after ensemble av- 
eraging. 

For odd n, T has no consequences because it has no 
eigenstates; since = —1, the (1 ±T)|JM) combina- 
tions change sign under T. 
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To (try to) explain JOD we note that all the J = 
states of seniority zero are in the self-conjugate space. It 
is fairly simple to calculate their centroids and widths. If 
they contribute substantially to the totals, much of the 
J-OES will be due to them and they will concentrate on 
the overall J = averages thus explaining their special 
status. A strong hint in favor of this argument is the 
massive presence of seniority zero states indicated by the 
persistence of pair-transfer coherence in J = ground 
states Q. To sum up all the indications related to T: 

Whatever separates even and odd J is due to the self- 
conjugate space, and whatever is special about J = is 
due to the seniority zero states, entirely contained in this 
space. And, remember, T invariance does not imply JOD, 
it only suggests its frequent observation. 
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FIG. 1. (p/) centroids and widths. See text 

Now do 50 runs for n = 8, 9 in the pf shell, normaliz- 
ing each to unit scalar variance cr^Jn = 2) = 1 . (For the 
calculation of traces we refer to ||ll|-[l3|l, though here we 
use a home-made method.) 

In Fig. [| the dots show Ecj (below) and aj (above) 
for the 50 71 = 8 runs, while the — full, dashed and 
dotted — lines are the averages for n = 8,9, and the re- 
sults for the full KB3 interaction plus single particle field 
in ''^Ca — respectively. 

As expected, the E^j ensemble average is smooth and 
for a J the OES between spins (J-OES) for n = 8 persists 
after ensemble averaging while cr(9)j does not stagger; 
and though larger on the average than cr(8)j, its maxi- 
mum is below f7(8)j=o- 

The random runs were done with monopole-free inter- 
actions {J2j ^rsrsC^J + 1) = 0)7 which simplify consid- 
erably the calculation of aj with negligible loss of ac- 
curacy. The '*^Ca results are exact (all eigenvalues for 
all J were calculated |^). Note that — to within a trivial 
overall factor — KB3 with its carefully tuned monopoles 
and experimental single particle spacings yields the same 
pattern for aj as the monopole-free ensemble average. 
No assumptions on two body monopoles will be made in 
what follows. 



The next task is to prove that Ecj and aj are suffi- 
cient to determine the ensemble averaged yrast patterns. 
The proof is as follows: The low moments of H deter- 
mine a smooth tridiagonal matrix. Once diagonalized it 
leads to a smooth binomial that describes very well the 
level densities. The position of the yrast state depends 
on the parameter Nj — hidj/ln2. If the third moment 
vanishes — as assumed here — the energy converges to the 
lower bound Ebj = Ecj — y^Njaj For moderate 

values of Nj (« 10 in our example) this limit is missed 
by « 2a J / \J Nj and a good estimate demands explicit 
diagonalization, which will obviously reflect J-OES. 

This result is strictly true for ensemble averages. For 
individual matrices, because of fluctuations, the bound 
is only an unreliable estimate (an old story under a new 
guise [§|15|). For example: with KB3, in ''^Ca the J 7^ 
yrast states come below their exact positions by 0.25 to 
2.5 MeV, while the ground state comes 4 MeV too high. 

The subject demands a full treatement, but the ar- 
guments above are sufficient to confirm the crucial role 
of a J anticipated in ^ but ignored in and put in 
doubt in through a misunderstanding that deserves 
a comment: Boson simulations are conducted in a collec- 
tive subspace of much smaller dimensions than the cor- 
responding fermion problem. In particular, the Nj x Nj 
tridiagonals in jl^ (once corrected |l^) are interesting 
models for the Lanczos submatrix at the origin , but it 
makes no sense to compare their variance to that of the 
full dj X dj matrix {dj = 2^-'). 

Of the three general properties mentioned in the first 
paragraph, BE2 enhancement is the one not sponta- 
neously produced by purely random trials. Some co- 
herence is needed, and the only way to simulate it, in 
a GOE context, is through a constant displacement, c, 
of all matrix elements. The displaced GOE (DGOE) is a 
standard ensemble, for which there exists a famous result: 
its spectrum is a semicircle, with one detached level 
Therefore, for attractive forces, the coherent part is a 
matrix whose elements are constants c = — |c|; leading to 
a displaced TBRE (DTBRE). 

In a two-body context with good J, attractive forces 
do not have systematically negative matrix elements, but 
their signs must have a very general origin, because all re- 
alistic interactions are spectacularly similar [|l0| , and the 
extremely rare sign discrepancies only affect the smallest 
absolute values. Furthermore, as was noted in Ref. [ p^ , 
these signs are strongly correlated to those of Elliott's 
q ■ q force. Hence, there may be a — DTBRE vs q ■ q — 
"sign coherence" conflict. 

Quite conveniently, it is possible to analize — and then 
resolve — it, while staying within the DTBRE; by exam- 
ining results in two spaces. The first is a (Aj = 2) sub- 
space of a major shell, in which aW q ■ q and realistic 
matrix elements are negative. It consists of the orbits 
with j = 1 + 1/2, for which a quasi- S't/S symmetry can 
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operate, leading to quadrupole properties similar to those 
of the full shell [Q. The second is the major shell itself. 

Accordingly, two sets of runs were performed. The first 
in the Aj = 2, (/7/2P3/2) space {fp for short) respects 
sign coherence. The other, in the full pf shell, does not. 
Please do not confuse fp and pf. 

960 fp runs were done for each of four combinations 
of number of particles and isospin: nT — 84(''*Ca), 40 
(^^Ti), 61 (4«^Ti), and 80 C'^Cr) , with strict GOE interac- 
tions plus single particle splitting. The parameters were 
chosen to mock realistic values {v = 0.6, e/^j^^ ~^h/2 ~ ^ 
MeV). The c = -1 steps for DTBRE are arbitrary. 

Fig. H shows the evolution as a function of c of the 
R — Ei^jEi energy ratio, a time-honored indicator of 
collective behaviour: R — \, 2, 3.33 corresponding to 
seniority, vibrational and rotational regimes respectively. 
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FIG. 2. Probability density for energy ratios E4/E2 

Only runs with J — ground states were kept and the 
probability density P normalized to unity for the four 
nuclei. At c = 0, the peaks are centered at i? = 1 ex- 
cept for ^^Ti. As c decreases ^^Ca will not budge, while 
for its partners R = 3.33 will be either clearly or over- 
whelmingly favored. There is more than a hint of rota- 
tional behaviour here, as confirmed by the distributions 
of BE2{2 0) strength in Fig. |. It is seen that at 
c = there is no coherence: the probabilities are con- 
sistent with a Porter-Thomas law. At c = — 1 the rates 
increase sharply, and at c = —2,-3 they concentrate 
in narrow peaks, whose strength should be compared to 
that of a rotor of intrinsic quadrupole moment Qo, and 
BE2{2 ^ 0) = 0.02^2/3 Ql e^Jm^. 

Following Ref. j2^ , we calculate the maximum Qo val- 
ues for 44Ti, 46Ti, and ^^Cr in the fp space to be (22, 25, 
32) (a-dimensional oscilator values) leading to transitions 
of (120, 160, 270) e^/m'* respectively. The main peaks 
are seen to be very close to these values, confirming the 
rotational nature of the T = 0, and 1 spectra. 

To check the influence of sign coherence, some 200 runs 
for "^^Cr were performed in the full pf shell, with uni- 
form single particle spacings of 2 MeV (as seen in ^^Ca). 



The similarities with the sign-coherent fp trends are as 
significant as are the differences. In both cases, as c be- 
comes more negative there is a consistent buildup of BE2 
strength, and a consistent increase of J = ground states 
(see Table | for the fp runs), increasingly associated to 
perfect yrast sequences (i.e, J=0,2,4,6,...). 
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FIG. 3. Probability densities for BE2{2 0) 

The discrepancies can be summed up by saying that in 
the full pf shell nothing is as clearcut as in the fp space, 
especially the limiting behavior for large displacements. 

The situation at c = — 2 is typical. In the pf runs we 
find "only" 75% of J = ground states (against 99.8% 
in fp), and practically as many perfect sequences. The 
indicators are good — but more spread than in Figs. ^ 
and |— at E4/E2 = 2.92 ± 0.56 and B£:2(2 ^ 0) = 
276±23 e^/m'' (the SU3 limit is 320 e^/m^, the strength 
is equally large in J 7^ ground states). More dis- 
turbingly, the runs do not seem to converge to perfect 
rotors: c = —3 brings little change over c — —2. The 
fairly good rotational behaviour suggested by the aver- 
ages above is not systematically confirmed by a well de- 
fined intrinsic state. In other words: Qo is not always 
approximately constant for the lowest members of the 
band. In the fp case, at c = —2, 99.8% of the sequences 
are perfect and all yrast levels are those of a rotor. 

Statistically, sign-coherence is not indispensable to 
generate acceptable BE2 enhancements, but physically 
it matters: Even nuclei do not have J = ground states 
50 or 80% of the time, but 100% of the time. 
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TABLE I. Percentage of J = ground states as a function 
of displacement c. 
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^^Cr has better than 75% chances of having a 
rotational-hke behaviour: it is a backbending rotor. 

Fig. m illustrates what a constant interaction W^'^f^ = 
—a does in the Aj = 2, {gds)^T = space. With single 
particle spacings of 1 MeV, the Ej+2 — Ej patterns are 
those of backbending rotors. They show more structure 
than the realistic interaction results but the physics 
is very much the same. 
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FIG. 4. Backbending patterns in {gdsf T = 0, with KLS 
interaction from and constant W^^^ ~ ~a- 

At this point, as far as DTBRE in the pf shell goes, 
we are where we we were twenty years ago in the sd shell: 
The BE2 coherence is there but "The spectra are neither 
rotational nor particularly interesting" | p9| . But now we 
have the fp spectra: rotational, and therefore interest- 
ing. The difference is obviously due to sign-coherence. 
But: Why do we have sign-coherence in fp7 Because in 
Aj = 2 spaces, the signs, by construction, are the same 
as in LS coupling j2^. 

The heart of the problem is that we do, routinely, 
simulations in jj coupling because it is the one used in 
shell model codes, as commanded by the central — more 
generally monopole — field. However, the "residual" two 
body multipole part of the Hamiltonian is dominated by 
forces that are overwhelmingly central (quadrupole, pair- 
ing, etc.) For these, the natural coupling scheme is 
LS. Therefore, it is artificial to define for them DTBRE 
coherence through the jj matrix elements. 

The solution of the sign-coherence problem becomes 
evident: define the DTBRE in LS scheme. Technically, 
we can continue to employ jj scheme using the LS dis- 
placement calculated once and for all. It is a safe bet — 
backed by the sd experience — that the pf results will 
become as satisfactory as the fp ones. 

Our examples involve rotational motion, but now it 
should be possible to generate other forms of collectiv- 
ity, as seen by the simplest example: with the same two 
body interaction the (sd)'* spectrum is rotational (^'^Ne) 
and the (sd)"'* spectrum is vibrational ('^^Ar). The dif- 
ference between the two is entirely due to the change of 



central field. Therefore, to obtain both rotors and vibra- 
tors it seems convenient to randomize the single particle 
energies and fix the two body terms, contrary to what is 
usually done. 

We have identified time reversal invariance as the ori- 
gin of odd-even staggering of the mass surfaces and J = 
spin for the ground states of even nuclei. The attractive 
nature of the forces appears to provide a sufficient condi- 
tion BE2{2 0) enhancements associated to collective 
behaviour. 
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